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Abstract 

In several extensions of the Standard Model, the top quark can decay into a bottom quark and a 
light charged Higgs boson H + , t -> bH + , in addition to the Standard Model decay t -> 6W. Since 
W bosons decay to the three lepton generations equally while H + may predominantly decay into 
tv, charged Higgs bosons can be searched for using the violation of lepton universality in top quark 
decays. The analysis in this paper is based on 4.6 fb _1 of proton-proton collision data at ^ = 7 TeV 
collected by the ATLAS experiment at the Large Hadron Collider. Signatures containing leptons (e or 
/i) and/or a hadronically decaying r (r had ) are used. Event yield ratios between e+r had and e+fi, as well 
as between ^ + r ha d and ^ + e, final states are measured in the data and compared to predictions from 
simulations. This ratio-based method reduces the impact of systematic uncertainties in the analysis. 
No significant deviation from the Standard Model predictions is observed. With the assumption that 
the branching fraction B{H + -»■ tv) is 100%, upper limits in the range 3.2%-AA% can be placed on 
the branching fraction B(t ->• bH + ) for charged Higgs boson masses m H + in the range 90-140 GeV. 
After combination with results from a search for charged Higgs bosons in it decays using the T had +jets 
final state, upper limits on B(t ->■ bH + ) can be set in the range 0.8%-3.4%, for m H + in the range 
90-160 GeV. 
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1 Introduction 



Several non-minimal Higgs scenarios, e.g. Two Higgs Doublet Models (2HDM) [1] predict 
the existence of charged Higgs bosons (H + and if"). 1 Their observation would clearly indi- 
cate physics beyond the Standard Model (SM), because this theory contains no elementary 
charged scalar particle. In several models, e.g. a type-II 2HDM describing the Higgs sector 
of the Minimal Supersymmetric extension of the Standard Model (MSSM) [2-6], the main 
production mode for charged Higgs bosons with a mass mjj+ smaller than the top quark 
mass (mtop) is through top quark decays t — > bH + . The dominant source of top quarks at 
the Large Hadron Collider (LHC) is through ti production. 2 

The combined LEP lower limit on the charged Higgs boson mass is about 90 GeV [7]. 
Results from direct searches for charged Higgs bosons decaying via H + — > tv using 4.6 fb _1 
of LHC data were recently presented by the ATLAS collaboration [8], with upper limits on 
the branching fraction B{t — > bH + ) between 5% and 1% for charged Higgs boson masses 
ranging from 90 GeV to 160 GeV, respectively. Using about 2 fb _1 of LHC data, the CMS 
collaboration established upper limits on B(t -> bH + ) in the range 4-2% for charged Higgs 
boson masses between 80 GeV to 160 GeV [9]. In all of these measurements, as well as in 
this paper (unless otherwise specified), the assumption B(H + — > tv) = 100% is made. This 
paper uses an alternative technique for H + searches in the mass range 90-160 GeV. In top 
quark decays, W bosons decay equally to leptons of the three generations, while H + may 
decay predominantly into tv. Hence, an excess of ti events with at least one hadronically 
decaying r lepton (rh a d) in the final state, as compared to the rate for ti events with only 
electrons and/or muons, is a signature for charged Higgs bosons. In order to cancel most 
of the systematic uncertainties in the analysis, a measurement of event yield ratios Ri for 
ti —> bb + /rhad + Nv and ti —> bb + IV + Nv, where Nv stands for any number of neutrinos 
and where I and V are electrons and muons, with I ^ V, is performed: 

B(ti^ bb + W hSiA + Nv) 
1 B(ti -»• bb + IV + Nv) ' ( ' ' 

This study is performed in a model-independent way, and so exclusion limits are given 
in terms of B{t — > bH + ), as well as in the m™ ax scenario [10] of the MSSM. The results 
are based on 4.6 fb _1 of data from pp collisions at ^fs = 7 TeV, collected in 2011 with the 
ATLAS experiment [11] at the LHC. These data, as well as the simulated samples used 
in the analysis, are the same as in ref. [8] and are described briefly in section 2. Then, 
in section 3, an event selection aimed at collecting a data sample enriched in ti events 
is presented, together with the data-driven methods to estimate the backgrounds due to 
misidentified electrons, muons and hadronically decaying r leptons. Exclusion limits in 
terms of B{t — > bH + ) and tan/3 are discussed in section 4, based on the measured ratios 
among event yields in Thad+lepton and dilepton final states. Finally, a summary is given 
in section 5. 

1 In the following, charged Higgs bosons are denoted H + , with the charge-conjugate H~ always implied. 
Hence, r denotes a positively charged r lepton. 

2 Since the cross section for H + production from events containing a single top quark is much smaller, 
this production mode is not considered here. 
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2 ATLAS data and simulated events 

The ATLAS detector [11] consists of an inner tracking detector with coverage in pseudora- 
pidity 3 up to \rj\ = 2.5, surrounded by a thin 2 T superconducting solenoid, a calorimeter 
system extending up to \rj\ = 4.9 for the detection of electrons, photons and hadronic jets, 
and a large muon spectrometer extending up to \r/\ = 2.7 that measures the deflection of 
muon tracks in the field of three superconducting toroid magnets. A three-level trigger 
system is used, which reduces the recorded event rate to about 300 Hz. 

In ATLAS, electrons are reconstructed by matching clustered energy deposits in the 
electromagnetic calorimeter to tracks reconstructed in the inner detector, and muons are 
required to contain matching inner detector and muon spectrometer tracks. The combina- 
tion of all sub-systems provides precise lepton measurements in the pseudorapidity range 
\r/\ < 2.5. Jets, as well as the magnitude of the missing transverse momentum, i?™ lss , 
are reconstructed using energy deposits over the full coverage of the calorimeters, up to 
\rj\ = 4.9. The anti-/cj algorithm [12, 13] with a radius parameter value of R = 0.4 is used for 
jet reconstruction. At least 75% of the tracks associated to a jet (weighted by their trans- 
verse momenta) must point to the primary vertex. This requirement on the "Jet Vertex 
Fraction" [14] allows the identification of jets originating from the hard-scatter interaction. 
An algorithm combining impact-parameter information with the explicit observation of a 
secondary vertex [15] is used in order to identify jets initiated by 6-quarks, within \rj\ < 2.4. 
The working point chosen for this study corresponds to an average efficiency of about 70% 
for 6-jets with > 20 GeV in tt events and a rejection factor of about 130 for light-quark 
jets. In order to reconstruct hadronically decaying r leptons, anti-fct jets with either one 
or three associated tracks, depositing Et > 10 GeV in the calorimeter, are considered as 
r candidates [16]. Dedicated algorithms are used to reject electrons and muons. The r 
candidates are further required to have a visible transverse momentum pL > 20 GeV and 
to be within |r/| < 2.3. The hadronic r decays are identified using a likelihood criterion 
designed to discriminate against quark- and gluon-initiated jets. The working point chosen 
for this study corresponds to an efficiency of about 30% for hadronically decaying r leptons 
with > 20 GeV in Z — > tt events, leading to a rejection factor of about 100-1000 for 
jets. Selected r candidates fulfilling the identification criteria are referred to as "t jets". 

The same full 2011 data set and simulated samples as the analysis in ref. [8] are used, 
corresponding to an integrated luminosity of 4.6 fb , with an uncertainty of 3.9% [17, 18]. 
In addition to the SM pair production and decay of top quarks, tt — > bbW + W~ , 4 the 
background processes include the production of single top quark, VF+jets, Z/7*+jets, 
diboson, and multi-jet events. Except for the last, which is estimated using data-driven 
methods, the SM backgrounds are determined using the simulated samples summarised 
in table 1. In addition to the SM background samples, three types of signal samples 

3 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The z-axis points from the IP to the centre 
of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, (j>) are used in the transverse 
plane, <f> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 9 as rj = — lntan(#/2). 

4 In the simulated SM ti events, all leptonic W decay modes have the same branching fraction (10.8%). 
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are produced with PYTHIA 6.425 [19] for 90 GeV < m H + < 160 GeV: tt -> bbH + W~, 
ti —> bbH~W + and ti — > bbH + H~, where charged Higgs bosons decay as H + — > tv. 
TAUOLA 1.20 [20] is used for r decays, and PHOTOS 2.15 [21] is used for photon radiation 
from charged leptons. All generated events are propagated through a detailed GEANT4 
simulation [22, 23], and they are reconstructed using the same algorithms as the data. 



Process 


Generator 


Cross section [pb] 


SM ti with at least one lepton I = e, fx, r 


MC@NLO 4.01 [24] 


91 


[25] 


Single top quark i-channel (with £) 


AcerMC 3.8 [26] 


21 


[27] 


Single top quark s-channel (with £) 


MC@NLO 4.01 [24] 


1.5 


[28] 


Single top quark VTt-channel (inclusive) 


MC@NLO 4.01 [24] 


16 


[29] 


W -+lv 


ALPGEN 2.13 [30] 


3.1 x 10 4 


[31] 


Z/j* -> ££ with m(££) > 10 GeV 


ALPGEN 2.13 [30] 


1.5 x 10 4 


[32] 


WW 


HERWIG 6.520 [33] 


17 


[34] 


zz 


HERWIG 6.520 [33] 


1.3 


[34] 


wz 


HERWIG 6.520 [33] 


5.5 


[34] 


H + signal with B(t -> bH + ) = 3% 


PYTHIA 6.425 [19] 


9.9 





Table 1. Cross sections for the simulated processes and the generators used to model them. 
All background cross sections are normalised to next-to-next-to-leading-order (NNLO) predictions, 
except for diboson event production where the next-to-leading-ordcr prediction (NLO) is used. For 
the diboson events, a filter is applied at the generator level, by requiring at least one electron or 
muon with px > 10 GeV and \r}\ < 2.8. 



3 Event selection and background determination 

This analysis uses events passing a single-lepton trigger with an Et threshold of 20 GeV or 
22 GeV for electrons and a px threshold of 18 GeV for muons. In order to select a sample 
enriched in ti events, the following requirements are made: 

• one charged lepton I (e, fj) having Er > 25 GeV (e) or px > 25 GeV (/x) and matched 
to the corresponding trigger object; 

• at least two jets having px > 20 GeV and \rj\ < 2.4, including exactly two 6-tags; 

• either exactly one r jet with p T > 25 GeV and \rj\ < 2.3 with no additional charged 
lepton, or exactly one additional charged lepton I' with £x or px above 25 GeV and 
a different flavour than the trigger-matched lepton; 

• E™ iss > 40 GeV. 

At this stage, the selected events are classified into two categories according to the 
single-lepton trigger that they fire: an electron trigger (EL) or a muon trigger (MU). Each 
category contains Th a d+lepton and dilepton (IV) events. The lepton appearing first in the 
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final-state is, by convention, matched to the corresponding trigger object. The EL category 
therefore consists of e + Th a d and e + fi events, while the MU category contains fi + r^ad 
and fi + e events. Events firing both a single-electron trigger and a single-muon trigger are 
assigned to both categories, and accounted for in the combined limit setting, see section 4.3. 

The analysis uses the generalised transverse mass mt^ 2 [35] as a selection variable. By 
construction, it gives an event-by-event lower bound on the mass of the charged (W or 
Higgs) boson produced in the top quark decay. Hence, it is larger than the true charged 
Higgs boson mass and smaller than m t0 p- For incorrect pairings of r jets or leptons 

with 6-jets, the numerical determination of m^ 2 may fail, hence only events with m^ 2 > 
are kept in the following. 

3.1 Backgrounds due to misidentified electrons and muons 

A significant background for the search described in this paper consists of events with 
reconstructed electrons and muons arising from the semileptonic decay of hadrons with b- 
or c-quarks, from the decay-in-flight of tt or K mesons and, in the case of electrons, from 7r° 
mesons, photon conversions or shower fluctuations. These are referred to as "misidentified 
leptons" in the following. Two data samples are defined, which differ only in the lepton 
identification criteria. The tight sample contains mostly events with real leptons and uses 
the same lepton selection as in the analysis. The loose sample is obtained by loosening the 
isolation and identification requirements, and it contains mostly events with misidentified 
leptons. 5 The efficiencies p r and p m for a real or misidentified lepton, respectively, to be 
detected as a tight lepton, are determined from data, with the same method as in ref. [8]. In 
the final parameterisation of p r and p m , dependencies on the pseudorapidity of the lepton, 
its distance AR = \J Ar] 2 + Acj) 2 to the nearest jet and the leading jet px are taken into 
account. Based on these efficiencies, the number of misidentified leptons passing the final 
requirements can be calculated by weighting each event in the data sample with one loose 
lepton, according to the following per-lepton weights wf. 

PmPr 



for a loose but not tight lepton, wn, 
for a tight lepton, wn 



Or - Pm) ' 

PmiPr ~ 1) 



(Pr ~ Pm) 

3.2 Backgrounds due to misidentified r jets 

About 51% of the simulated ti events in the Th a d+lepton final state contain a r jet matched 
to a hadronically decaying r lepton at the generator level. In the other events, the r jet 
is called "misidentified" . It originates from leptons (e, fj,) in 3% of the simulated events 
and hadronic objects (initiated by light quarks, 6-quarks or gluons) in 46%. Data-driven 
methods are used in order to determine the probability of misidentification in the case 
of hadronic jets and electrons. In the latter case, the misidentification probabilities are 
measured using a Z — > ee control region in the data [16] and then applied to the simulated 
events, as in the analysis in ref. [8]. 



3 By construction, the tight sample is a subset of the loose sample. 
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The majority of misidentified r jets in the final event selection originate from jets, for 
which the misidentification probability depends on the initial parton (light quark, heavy- 
flavour quark or gluon). All jet types occur in ti events, and it is not possible to accurately 
predict the fraction of each of them, potentially leading to a large systematic uncertainty 
on the jet — > Th a( j misidentification probability. However, the influence of all jet types 
other than light-quark jets can effectively be eliminated by categorising all events in terms 
of the charge of the lepton relative to the r jet as opposite-sign (OS) or same-sign (SS) 
events. All processes with gluon and 6-quark jets produce positively and negatively charged 
misidentified r objects at the same rate. On the other hand, the light-quark jet component 
in SS events represents both charge misreconstruction and quarks which fragment such that 
the leading charged particle does not have the same charge as the initial quark. Giving a 
negative weight to the SS events therefore cancels, on average, the gluon and heavy-flavour- 
quark jet contributions from the OS events, leaving only light-quark jets misidentified as 
r jets. 

The rate at which light-quark jets are misidentified as r candidates is derived using a 
region enriched with W + >2 jets events in the data, selected by requiring: 

• exactly one electron or muon with Ej> or larger than 25 GeV; 

• at least one r candidate; 

• at least two jets in addition to the r candidate(s), none of them being 6-tagged; 

• Elf iss > 40 GeV. 

The leading process in this control region is gq —> Wq'. Hence, the final state quark 
misidentified as a r candidate and the W boson usually have opposite charges. In order 
to reduce the contribution from events with a true r lepton, mostly from Z+jets events, a 
requirement on the transverse mass wt is made: 



where A(^ m i ss is the azimuthal angle between the lepton and the direction of the missing 
momentum. The W + >2 jets events are classified as OS and SS events using the charges 
of the lepton and the r candidate. Figure 1 shows the wit distribution for OS, SS and 
OS-SS events fulfilling the W + >2 jets selection. This demonstrates the cancellation of 
heavy-flavour-quark and gluon contributions. 

However, figure 1 also suggests that simulation poorly models r candidates, which in 
turn motivates the use of data-driven methods to predict the background due to misiden- 
tified r jets. In particular, the number of tracks associated to jets misidentified as r 
candidates is found to be poorly modelled in simulation. Events in the data tend to have 
fewer r candidates with one or three tracks. In order to correct the r candidate selection 
efficiencies in simulation, r track multiplicity scale factors are derived using OS-SS events 
fulfilling the W + >2 jets selection, and are then applied to all jets misidentified as r 
candidates in the simulation: 0.71 ± 0.03 for 1-track r candidates; 0.92 ± 0.03 for 3-track 
r candidates, where the errors are only statistical. 
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Figure 1. Distributions of the transverse mass mx for events fulfilling the W + >2 jets selection, 
without the requirement tot > 30 GcV. Each colour corresponds to a type of generator-level particle 
matched to a r candidate (i.e. the highest energy particle within a cone of radius Ai? = 0.2 around 
the r candidate is considered) . The SS events are (a) given a weight of — 1 and (b) subtracted from 
the OS events. All simulated SM processes are considered. 

The probability for a light-quark jet to be misidentified as a r jet is measured in the 
data and is binned in pT^, the number of associated tracks A^ ack (one or three), and the 
number of tracks Nl^° k found within 0.2 < Ai? < 0.4 of the r candidate. For each bin, 
the jet — > Thad misidentification probability is defined as the number of objects passing the 
r identification based on the likelihood criterion divided by the number prior to requiring 
identification. OS events are given a weight +1 and SS events are given a weight —1, in 
both the numerator and denominator of the jet — > Thad misidentification probability. After 
OS-SS subtraction, the selected events mostly contain r candidates coming from light- 
quark jets and, to a much lesser extent, electrons, muons, and true hadronically decaying r 
leptons. Figure 2 shows the measured values of the jet—)- Thad misidentification probability 
in W + >2 jets events selected from the data, after OS-SS subtraction. These are used to 
scale all simulated events in the signal region. Events fulfilling the requirements listed in 
the beginning of this section, in which the selected r object originates from a jet (of any 
type), are weighted by the misidentification probabilities. An additional weighting factor 
(+1 for OS events and —1 for SS events) is then used to perform the OS-SS subtraction. 

4 Results 

4.1 Computation of event yield ratios 

For each of the four final states considered here (e + Th a d) e + /x, /x + Th a d an d fJ> + e), 
the OS-SS event yield M can be split into two contributions: from ti events (where the 
top quark decays are mediated by both W and H + ) and from all other SM processes 
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Figure 2. Probability for a light-quark jet to be misidentified as a 1-track or 3-track r jet, measured 
in a region enriched with OS-SS W + >2 jets events in the data, as a function of (a) and (b) 
the number of tracks N£° ck found within 0.2 < AR < 0.4 of the r jet. 

except ti — > bbW + W~ . The contributions from ft events are expressed as a function of the 
cross section a t T, the integrated luminosity C, the branching fraction B = B(t — > bH + ), 
as well as the selection efficiencies e w + w -, e H + w -, e H - w + and e H + H - for, respectively, 
ft -> ft -> bbH + W~, ti -> bbH~W + and ft -> bbH + H~ events, in each of the 

four final states considered here: 

A/" = <7 tt - x £ x [(1 - B) 2 e^+W- + 5(1 - 5) (e H+w - + e a -w+) + £ 2 eff+#-] 

+ AWers • (4.1) 

In turn, event yield ratios are defined as: 

R e = ^±^4 ^d R» = ^ + Th ^ . (4.2) 

The event yields in the Thad+repton and dilepton final states are summarised in table 2 
for the background-only hypothesis, as well as in the presence of a 130 GeV charged Higgs 
boson in the top quark decay. The predicted values in the SM-only hypothesis and the 
measured values of the ratios R e and i? M are summarised in table 3. 

Figure 3 shows the variation of the event yields M{e + Th a d)> M{e + /i), N{n + Th a d) 
and ftf(fi + e) with B(t — > bH + ), for a charged Higgs boson mass of 130 GeV. The presence 
of H + — > tv in a fraction of the top quark decays leads to an increase of the number of 
ti events with a lepton and a r jet, hence of the ratios R e and R^. Note that the event 
yields for dilepton final states become smaller in the presence of a charged Higgs boson in 
top quark decays, despite the fact that a r lepton decays into an electron or muon more 
often than a W boson. This results from the fact that electrons and muons produced in 
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Sample 


OS-SS event yields 
e + Thad e + (j, 


Misidentified electrons or muons 
W/Z+}ets & diboson 
Single top quark 
ti 


-0.8 ± 3.0 94 ± 37 
2.1 ± 0.9 0.7 ± 0.4 
3.3 ± 0.8 24 ± 4 

111 ±25 980 ±200 


£ SM 


116 ±25 1100 ±210 


Data 


144 1247 


tt with t -> bH+ (130 GeV) 
Prediction with signal 


30 ± 4 27 ± 4 
139 ±28 1070 ±200 
M + T had M + e 


Misidentified electrons or muons 
W/Z+jets k diboson 
Single top quark 
tt 


0.2 ± 1.0 74 ± 37 
2.6 ± 1.6 0.7 ± 0.4 
4.6 ± 0.9 18 ± 3 
131 ±28 740 ±150 


£ SM 


138 ±29 830 ±160 


Data 


153 929 


it with t -> (130 GeV) 
Prediction with signal 


35 ± 4 20 ± 3 
166 ±32 810 ±150 



Table 2. Expected OS-SS event yields after all selection cuts in Thad+lepton and dilepton channels, 
compared with 4.6 fb _1 of ATLAS data. The numbers shown for a hypothetical 130 GeV H + signal 
correspond to B(t — > bH + ) = 3%. The contribution of tt — > bbWW events to the background is 
scaled accordingly Statistical and systematic uncertainties are combined. 



Ratio 


Re 


Rfi 


SM value 


0.105 ±0.012 


0.166 ±0.017 


Measured value 


0.115 ±0.010 (stat) 


0.165 ±0.015 (stat) 



Table 3. Predicted (in the SM-only hypothesis) and measured values of the event yield ratios 
R e and i? M . For the values of the ratios predicted using simulation, the statistical and systematic 
uncertainties are combined. 

the decay chain t — > bH + — > bru — )• bl + Nu are, on average, softer than those coming from 
t — > bW — > bl + Nv. The sensitivity of this analysis to charged Higgs bosons is determined 
by the rate at which the ratios R e and R^ change with B(t — > bH + ), which depends on 
the selection efficiencies e H + w -, e H - w +, e H + H - and, in turn, on the charged Higgs boson 
mass. For m H + = 150 (160) GeV, the rate at which the ratios R e and R^ change with 
B(t — > bH + ) is found to be two (five) times smaller than for rrijj+ = 130 GeV. Indeed, the 
selection efficiencies e H + w ~, e H ~ w +, e H + H - are reduced for m H + values in the vicinity 
of mtop, because the 6-jet arising from t — > bH + becomes softer when the mass difference 
"itop — r nH+ is smaller. 



-9- 




Figure 3. Relative variation with Bit — > bH + ) of (a) the event yields M{e + ri ia d), M{e + /x) and 
their ratio, as well as (b) 7V(/i + Th a d), Af(fJ. + e) and their ratio, assuming the presence of a 130 GeV 
charged Higgs boson in ti events. 

4.2 Systematic uncertainties 

Systematic uncertainties arise from the simulation of the electron and muon triggers, from 
the reconstruction and identification efficiencies of the physics objects, as well as from the 
energy /momentum scale and resolution for these objects. In order to assess their impact, 
the selection cuts of this analysis are re-applied after shifting a particular parameter by its 
±1 standard deviation uncertainty, while other parameters are fixed. The largest instru- 
mental systematic uncertainties are for jets. In comparison, the systematic uncertainties 
arising from the reconstruction and identification of electrons, muons and r jets are small. 
All instrumental systematic uncertainties are propagated to the reconstructed E™ ms . 

The ti cross section used in this analysis is a t i = 167l}g pb [25]. To estimate systematic 
uncertainties arising from the ti generation and the parton shower model, the acceptance 
is computed for ti events produced with either MC@NLO interfaced to HERWIG and 
JIMMY [36] for the hadronisation and the underlying event, or POWHEG [37] interfaced to 
PYTHIA. Systematic uncertainties on initial- and final-state radiation are computed using 
tt samples generated with AcerMC interfaced to PYTHIA, where the relevant parameters 
in PYTHIA are varied in a range given by current experimental data [38]. These systematic 
uncertainties are dominated by the difference in modelling the numbers of tracks N£ rack 
and N™° ck in, respectively, the core and isolation regions of the jets misidentified as r 
candidates. The various simulated ti samples are reweighted so that the iV£ ack and the 
-^tmck distributions match 6 before the systematic uncertainties on the ti generation, the 
parton shower model, as well as initial- and final-state radiation, are evaluated. 

6 Both variables are reweigthed in a correlated way. 
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For the signal samples, which are generated with PYTHIA (i.e. without higher-order 
corrections), no alternative generator is available, hence the systematic uncertainty is set to 
the relative difference in acceptance between ti events generated with MC@NLO interfaced 
to HERWIG/ JIMMY and with AcerMC, which is also a leading-order generator, interfaced 
to PYTHIA. For the systematic uncertainty coming from initial- and final-state radiation, 
the same simulated samples as for the SM ti events are used. In the evaluation of the 
systematic uncertainties for the signal samples, only r jets matched to true hadronically 
decaying r leptons in the generated events are considered. 

For the backgrounds with misidentified leptons, the largest systematic uncertainties 
arise from the sample dependence: the misidentification probabilities are calculated in 
a control region dominated by gluon-initiated events, but later used in a data sample 
with a higher fraction of quark- initiated events. The total systematic uncertainty on the 
backgrounds with misidentified leptons is 38% for electron-triggered events and 49% for 
muon-triggered events. It corresponds to the relative variation of the number of events 
with exactly one trigger-matched lepton and two jets, after having considered all systematic 
uncertainties. The requirement of having two 6-jets in the event does not have a significant 
impact on these systematic uncertainties and neither does the presence of a second lepton. 

For the estimation of backgrounds with jets misidentified as hadronically decaying r 
leptons, the systematic uncertainty on the scale factors associated with the number of 
tracks is determined by varying the requirement on the jet multiplicity and the magnitude 
of the subtraction of r candidates matched to a true electron, muon or r lepton in the 
generated events. This uncertainty is 7% for 1-track r jets and 11% for 3-track r jets. 
In addition, systematic uncertainties on the jet—)- r^ad misidentification probability arise 
from statistical uncertainties due to the limited control sample size, the differences between 
misidentification probabilities computed in the region enriched with W + >2 jets events 
and the signal region, as well as the small contamination from true r leptons (including 
those possibly coming from H + — > tu) in the region enriched with W + >2 jets events. 

Some of the systematic uncertainties above affect the Th a d+l e Pt° n an d dilepton event 
yields in the same manner and, as a result, have a limited impact on R e and R^. Systematic 
uncertainties arising from jets and I?™ 1SS are common to all reconstructed events in the 
simulation, hence they should cancel in the ratios R e and R^. However, due to the use of 
data-driven background estimates and because of the removal of geometric overlaps between 
reconstructed objects, some of these systematic uncertainties still have a minor impact. In 
the EL (MU) category, the systematic uncertainties related to the trigger-matched electron 
(muon) are the same for the e + Th a d and e + fj, (fj, + r^ad and // + e) events, thereby not 
affecting the predicted value of the ratio R e (R^)- Those coming from the reconstructed 
muon (electron) only affect event yields in the denominator, and hence the ratio. Similarly, 
the systematic uncertainties coming from the r jets and their misidentification probabilities 
only affect the numerator of R e and i?^, hence they do have an impact on the analysis. 
This is also the case for systematic uncertainties on the backgrounds with misidentified 
leptons, which have a larger contribution in the dilepton events, i.e. on the denominator 
of R e and R^. Table 4 shows how these ratios (in the SM-only hypothesis) change when 
shifting a particular parameter by its ±1 standard deviation uncertainty. 
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Systematic uncertainty 






Integrated luminosity 


0.3% 


0.3% 


Electron trigger efficiency 


0.1% 


N/A 


Electron reco. and ID efficiencies 


0.2% 


1.9% 


Electron energy resolution 


0.1% 


<0.1% 


Electron energy scale 




0.1% 


0.3% 


Muon trigger efficiency 




N/A 


0.1% 


Muon reco. and ID efficiencies 


1.0% 


0.1% 


Muon momentum resolution 


<0.1% 


<0.1% 


Muon momentum scale 




0.1% 


<0.1% 


r ID efficiency 




3.9% 


3.9% 


r energy scale 




2.9% 


3.0% 


r mis-ID (data-driven): 


number of associated tracks 


2.1% 


2.1% 


r mis- ID (data-driven): 


true Thad contamination 


0.2% 


0.2% 


r mis- ID (data-driven): 


H + signal contamination 


0.6% 


0.6% 


r mis-ID (data-driven): 


event environment 


1.3% 


1.2% 


r mis-ID (data-driven): 


statistical uncertainties 


3.3% 


3.2% 


r mis-ID (data-driven): 


electron veto uncertainties 


0.6% 


0.3% 


6-tagging 




1.9% 


2.3% 


Jet vertex fraction 




0.1% 


0.4% 


Jet energy resolution 




0.4% 


<0.1% 


Jet energy scale 




0.7% 


0.5% 


Jet reconstruction efficiency 


0.1% 


0.4% 


irimiss 


0.3% 


0.1% 


tt: cross section 




0.7% 


0.6% 


tt: generator and parton shower 


5.7% 


4.4% 


ti: initial- and final-state radiation 


3.6% 


3.7% 


Backgrounds with misidentified leptons 


3.5% 


4.3% 


Total (added in quadrature) 


10.3% 


10.1% 



Table 4. Relative variation of the ratios R e and in the SM-only hypothesis after shifting a 
particular parameter by its ±1 standard deviation uncertainty. 

4.3 Exclusion limits 

To test the compatibility of the data with the background-only or the signal+background 
hypotheses, a profile likelihood ratio [39] is used with R e and as the discriminating 
variables. The systematic uncertainties are incorporated via nuisance parameters, and the 
one-sided profile likelihood ratio, q^, is used as a test statistic. No significant deviation from 
the SM prediction is observed in 4.6 fb" 1 of data. Exclusion limits are set on the branching 
fraction B(t — > bH + ) by rejecting the signal hypothesis at the 95% confidence level (CL) 
using the CL S procedure [40]. These limits are based on the asymptotic distribution of 
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the test statistic [39]. They are first set for electron-triggered and muon-triggered events 
separately (see figure 4), and then using a global event yield ratio R e +^ defined as: 

AA(e + r had ) + Af(p + r had ) 



R 



e+fJ. 



(4.3) 



N{e + Li) + N OK {Li + e) ' 

where A/or(/^ + c) is the event yield in the Li + e channel after removing the dilepton events 
that simultaneously fire a single-electron trigger and a single-muon trigger, as those already 
appear in jV(e + li). The fraction of dilepton events common to the li + e and e + fj, final 
states is about 42% in the data. Using this global event yield ratio, upper limits in the 
range 3.2%-4.4% can be placed on B(t —> bH + ) for charged Higgs boson masses in the 
range 90-140 GeV, as shown in table 5 and figure 5. 
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Figure 4. Upper limits on Bit — > bH + ) derived from the event yield ratios (a) R e and (b) as 
a function of the charged Higgs boson mass, obtained for an integrated luminosity of 4.6 fb _1 and 
with the assumption B{H + — > rv) = 1. The solid line in the figure is used to denote the observed 
95% CL upper limits, while the dashed line represents the expected exclusion limits. The green 
and yellow regions show the ler and 2a error bands. 



m H + (GeV) 


90 


100 


110 


120 


130 


140 


150 


160 


95% CL observed 
(expected) limit on 
B(t-> bH+) using 
the ratio R e +n 


3.3% 
(3.1%) 


3.6% 
(3.3%) 


3.2% 
(3.0%) 


3.4% 
(3.1%) 


3.6% 
(3.3%) 


4.4% 
(4.0%) 


7.3% 
(6.7%) 


18.3% 
(16.8%) 



Table 5. Observed (expected) 95% CL upper limits on B(t —> bH + ) derived from the event yield 
ratio Re+ii, as a function of the charged Higgs boson mass, obtained for an integrated luminosity 
of 4.6 fb _1 and with the assumption that B(H + — > tv) = 1. 



In a previously published search for charged Higgs bosons [8], based on the data col- 
lected in 2011 with ATLAS, upper limits on B(t — > bH + ) were derived using various 
distributions of discriminating variables in Th a d+jets, Th a( j+lepton and lepton+jets final 
states. The most sensitive channel was Th a d+jets, except for low values of m H +. A new 
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Figure 5. Upper limits on B(t — > bH + ) derived from the event yield ratio R e +^, as a function 
of the charged Higgs boson mass, obtained for an integrated luminosity of 4.6 fb -1 and with the 
assumption B(H + — > tv) = 1. The solid line in the figure is used to denote the observed 95% CL 
upper limits, while the dashed line represents the expected exclusion limits. The green and yellow 
regions show the la and 2a error bands. 

set of combined upper limits on B(t — > bH + ) is derived, using the transverse mass distri- 
bution of Thad+jets events from ref. [8] and the global event yield ratio R e +^, as shown in 
table 6 and figure 6. Since a lepton veto is applied for charged Higgs boson searches in 
Thad+jsts final states, there is no correlation between such events and those selected in this 
study to determine the event yield ratios. With this combination of upper limits, charged 
Higgs bosons can be excluded for values of the branching fraction B(t — > bH + ) larger than 
0.8% to 3.4%, for m H + between 90 GeV and 160 GeV. These exclusion limits represent an 
improvement with respect to those published in ref. [8]. 
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(expected) limit on 


3.4% 


2.9% 


1.7% 
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1.1% 


B(t^bH+) using 
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(1.2%) 


R e+ll and T had +jets 


















95% CL observed 


















(expected) limit 


4.8% 


3.4% 


2.1% 


1.3% 


1.1% 


1.0% 


1.1% 


1.0% 


in ref. [8] 


(4.2%) 


(3.5%) 


(2.5%) 


(1.9%) 


(1.5%) 


(1.3%) 


(1.2%) 


(1.3%) 



Table 6. Observed (expected) 95% CL upper limits on B(t — > bH + ) derived using Thad+jets 
events in ref. [8] and the ratio R e + l n as a function of the charged Higgs boson mass, obtained for 
an integrated luminosity of 4.6 fb -1 and assuming that B(H + — > tv) = 1. The exclusion limits 
published in ref. [8] are also shown for comparison purposes. 
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Figure 6. Upper limits on Bit — s- bH + ) derived from the transverse mass distribution of Th a d+jcts 
events in ref. [8] and the event yield ratio R e +^., as a function of the charged Higgs boson mass, 
obtained for an integrated luminosity of 4.6 fb -1 and with the assumption B(H + — > tv) = 1. The 
solid line in the figure is used to denote the observed 95% CL upper limits, while the dashed line 
represents the expected exclusion limits. The green and yellow regions show the la and 2a error 
bands. 

In figure 7, the combined limit on B(t — > bH + ) x B{H + — > tv) is interpreted in the 
context of the m™ ax scenario [10] of the MSSM. The following relative theoretical uncer- 
tainties on B(t —> bH + ) are considered [41, 42]: 5% for one-loop electroweak corrections 
missing from the calculations, 2% for missing two-loop QCD corrections, and about 1% 
(depending on tan /3) for A^-induced uncertainties, where is a correction factor for the 
running 6-quark mass [43]. These uncertainties are added linearly, as recommended by the 
LHC Higgs cross-section working group [42]. 

Assuming that the boson recently discovered at the LHC [44, 45] is one of the neutral 
MSSM Higgs bosons, only a certain region in the m^+-tan/3 plane is still allowed for a 
given scenario [46]. If the new boson is the lightest neutral MSSM Higgs boson 
it would imply tan/3 > 3 and rritf+ > 155 GeV. However, the allowed region depends 
strongly on MSSM parameters which, on the other hand, do not affect the charged Higgs 
boson production and decay significantly. Thus, by adjusting these MSSM parameters, 
the region in which the Higgs boson mass can take a value of about 125 GeV can be 
changed significantly, while the ATLAS exclusion region shown here is relatively stable 
with respect to these changes. Should the recently discovered boson instead be the heavier 
CP-even Higgs boson (H°), the additional constraint from m^o — 125 GeV only leads 
to an upper limit of roughly m H + < 150 GeV, with suppressed couplings for h°. If the 
recently discovered particle is an MSSM Higgs boson, excluding a low-mass charged Higgs 
boson would thus imply that it is the lightest neutral state h°. 
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Figure 7. Limits for charged Higgs boson production from top quark decays in the m^+-tan/3 
plane, derived using Thad+jets events in ref. [8] and the ratio i? e +/j, in the context of the m™ ax sce- 
nario of the MSSM. The la band around the observed limit (dashed lines) shows the theoretical un- 
certainties. Values below tan/3 = 1, where the calculations in the MSSM become non-perturbative, 
are not considered, as the results become unphysical. 

5 Conclusions 

Charged Higgs bosons have been searched for in it events, in the decay mode t — > bH + 
followed by H + — > tv. A total of 4.6 fb" 1 of pp collision data at -y/s = 7 TeV, recorded 
in 2011 with the ATLAS experiment at the LHC, is used. Event yield ratios are measured 
in the data and compared to the predictions from simulations, between electron-triggered 
e + Thad arid e + /i events, and between muon-triggered /i + Thad and /i + e events, in 
order to search for a violation of lepton universality in ti events. This method reduces 
the impact of several systematic uncertainties in the analysis. Data-driven methods and 
simulation are employed to estimate the number of background events. The observed data 
are found to be in agreement with the SM predictions. Assuming B(H + — > tv) = 100%, 
upper limits at the 95% confidence level in the range 3.2%-4.4% have been placed on the 
branching fraction B{t — > bH + ) for charged Higgs boson masses in the range 90-140 GeV. 
For charged Higgs boson masses below 110 GeV, this analysis improves the previously 
published limits on B{t — > bH + ), based on direct searches for charged Higgs bosons in ti 
decays using the lepton+jets, Thad+jets and Th a d+lepton final states. When the results 
of the present analysis are combined with the results from the search for charged Higgs 
bosons in ti decays using the Thad+jets final state [8], upper limits on B{t —> bH + ) are set 
in the range 0.8%-3.4%, for m#+ between 90 GeV and 160 GeV. 
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